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Abstract. We consider a discrete time particle model for kinetic transport on the two 
dimensional integer lattice. The particle can move due to advection in the x-direction and 
due to dispersion. This happens when the particle is free, but it can also be adsorbed and 
then does not move. When the dispersion of the particle is modeled by simple random 
walk, strange phenomena occur. In the second half of the paper, we resolve these problems 
and give expressions for the shape of the plume consisting of many particles. 
Key words, simple random walk, conditional variance, reactive transport, kinetic adsorp- 
tion, double-peak. 
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1. Introduction 

We consider a discrete time two dimensional stochastic kinetic transport model, which 
is described by a two dimensional random walk with a drift in the x direction. Here the 
kinetics refers to the possibility that the particle may be adsorbed and released at random 
times, slowing down the movement of the particle. Let S{n) — (Sxin), S'y(n)) be the position 
of the particle at time n — 0,1, ... . In the engineering literature (see e.g. simulations are 
performed with a large number of particles independently performing such a slowed down 
random walk. This results in a plume of particles, and there is interest in the shape of this 
plume. One of the remarkable phenomena is that level sets of the concentration function of 
the plume need not be convex (see Figure [l] last column). 

A natural quantitative way to measure the lateral spread of the plume at position x at 
time n is to consider cr^ — Var (S'y(n) | Sxin) ~ x). The goal of our paper is to compute 
this conditional variance, and to discuss its properties. Intuitively, cr^ should be increasing 
in X. However, if we model the random walk with the natural choice: the simple symmetric 
random walk on the 2D integer lattice, then (see Figure [l] bottom row) this is not the case, 
as we prove in Section [3] A second strange phenomenon that we observe here is that there is 
symmetry (i.e., cr^ = (t|„_j.) when the total number of times the particle has been free follows 
a binomial distribution (see Figure [l] bottom left). The background for this is a symmetry 
property of the asymmetric 2D simple random walk, as we show in Section [4j 
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In Section [5] we get rid of the strange behaviour caused by the simple random walk by 
tilting the steps. We also show that this choice is 'right', as we obtain the same behaviour in 
the case of Gaussian distributions for the random walk. 

For related work lying at the basis of the present paper, see [S]. 




Figure 1. The probability mass function /50 of K^q, the probability mass 
function of 5(50), a list contour plot of its mass function and the conditional 
variance Var (5*^(50) | S'x(50) = x) with = tt, a = /3 = 1/4. In the first 
column a = 0.1, b = 0.9; in the second column a ~ h ~ 0.1; in the last 
column a = 6 = 0.01. See Section[2]and Example |3.1| for further explanations. 



2. A DISCRETE TIME STOCHASTIC MODEL FOR KINETIC TRANSPORT 

The kinetic transport model is built out of two independent processes: a 2D random walk 
{{Xk,Yk), k > 1}, and a two-state Markov chain {Zk, A; > 1} (see [2] for the ID case). We 
denote by P the probability measure on a large enough space to accomodate both processes. 
The two processes describe the behavior of a single particle that in each unit time interval 
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can only be in one of two states: 'free' or 'adsorbed', coded by the letters F and A. The 
particle can only move when it is 'free', and in this case its displacement has two components: 
dispersion and advection. 

Let {Xk,Yk) be the displacement of the particle due to the dispersion the fcth time that 
it is 'free'. We model the {Xk,Yk) as independent identically distributed random vectors 
satisfying 

(1) P((Xfe,yfc) = (j,0))-a, P((Xfe,n) = (0,j))=^ forj=±l 

where a, /? > 0, a + /3 — 1/2. This means that the particle moves one unit to the left or the 
right with the same probability a, and moves one unit up or down with the same probability 
/3. When the particle is free, the displacement during a unit time interval due to advection 
is given by one unit in the x direction. 

The kinetics, i.e., the movement from 'free' to 'adsorbed' and vice versa, is modeled by a 
Markov chain {Zk, /c > 1} on the two states {F, A} with initial distribution v = {vp, vx) and 
transition matrix 

1 — a a 
b 1-b ' 

Clearly the stationary distribution tt — {t:y,t^a) of the Markov chain {Zi^} is given by 7rp = 
b/{a + b) and tta — a/{a + b). Let Kn be the occupation time of the process {Zk} in state F 
up to time n, i.e., 

n 
k=l 

Then the distribution of Kn is well known, and is called a Markov binomial distribution (see, 
e.g., lU HI). Let /„ be its distribution, i.e., /n(fc) = P [Kn = k) ior < k < n. Then it is 
easy to calculate /„(fc) by the following recurrence equation 

/n+2(fc + !) = (!- b)U+iik + !) + (!- a)/„+i(fc) - (1 - a - 6)/„(fc) 
with initial conditions 

/i(0) = i^A; /i(l) - i^F, /2(0) = i^Ail - b), hil) = i^Ab + y^a, /2(2) - Ml - a). 

Now let S{n) be the position of the particle at time n, where 5'(0) = (0, 0). Then by the 
above we can write S(n) for n > 1 as 

S{n) = {Sx{n),SY{n)) = Y.^Xk + l,Yk). 

fc=i 

3. Conditional variance 

In this section we are interested in the conditional variance of Syin) given Sx{n) = x. 
For a real number c € M, let [cj be the largest integer less than or equal to c, and [c] be the 
least integer greater than or equal to c. For two real numbers c and d, let c V d := max{c, d} 
and c A d := min{c, d}. 



(2) 



Pf,f Pf,a 
Pa,f Pa.a 
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Proposition 3.1. For n>\, and < x < 2n, —n < y < n, we have for x = y mod 2 

F{S{n) = ix,y))= /„(fc) _ _ _ , 

fc=[a;/2] j=ow{x-k)-^ " ; vv I «;/ vv yj/ J J 

and P {S{n) = (x, y)) = for x ^ y mod 2. 

Proof. For < a; < 2n and — n < y < n, we consider the probability of the particle being at 
position (x, y) at time n, conditioned on — k. Let x {x, y, y) be the number of rightwards 
(leftwards, upwards, downwards) steps of the random walk 5*^ = [Xi + ■ ■ ■ + Xk + k,Yi + 
■ ■ ■ + Yfe). Suppose the walk makes x = j steps to the right. On the event {Sk — {x,y)} we 
then have 

^ x^j ( x^j 

X — X = X — k \x = j + k — X 

y-y^y | y ^ {x + y)/2-' j 

x + x + y + y^k [ y = {x - y)/2 - j. 

Necessarily, x = y mod 2. Since x, x, y, y > 0, it immediately follows that 

(x + yx — y\ 
max{0,x- k} <j < mm|^— , j. 

By a standard combinatorial analysis, and using Equation ([l]) we see that 
P {S{n) = [x, y)\K^^ k) = P (5,- = (x, y)) 

k\ a^3+k-xpx~2j 



E 



Clearly x < 2k, yielding that k > \x/2']. Multiplying by fn{k), and adding over k ranging 
from \x/2~\ to n, completes the proof of the proposition. □ 

Proposition 3.2. For n > 1 and < x < 2n, the conditional variance is given by 

lx/2\ 

E fnik) E (x-2j)0.,.+,_t,_,,>^^+'^-^(2/3)-2^- 

■\r f \ I Q f \ \ k={x/2-] j=OW{x-k) 

Var (S'y(n) Sxin) = x) = — . 

n lx/2} 
k=[x/2] j=OV{x-k) 

where „^f_ J = e<. s\ (m-e-sy. denotes a trinomial coefficient. 
Proof. In a similar way as in Proposition |3.H we obtain that 

n 

P(^xH=x) = ^/„(A:)P(5xW=a; I = k) 

k=0 

lx/2\ 

- E /"(^) E .u^k-^x.u,^2jy -''''-m-''- 

k=\x/2'] j=oy(x-k)-^' ''^ 
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By symmetry E[S'y(7i) | Sx{n) — x] = 0. So 

Var {Sy (n) \ Sx (n) = x) = ^ y^P {Syin) ^ y \ Sx (n) = x) , 



H = x mod 2 

—n<y<n 



and the proposition follows by using Proposition 3.1 in the following calculation 
J2 y^P{S{n)^ix,y)) 



= x mod 2 

-n<.y<.n 



—n<y<n 

n [_x/2\ x-2] , . 

= y fn{k) y -T7 ^Att ^a''+'^-^f^^-''y{x-2j-2s)^'' 

^ ^ ' ^ + k-x)\{x-2j)\ ^ ■' \ s 

" La;/2J , 

k={x/2'] j=OW{x~k) J \J < J \ J) 

where the first equality follows from the usual convention (™) = for ^ < or £ > m. □ 

Example 3.1. Let n = 50, a — (3 — 1/4, and the initial distribution v = t:. See Figure^ 

Note that the last column of Figure [T] shows that two peaks appear in the probability mass 
function of S{n) when both a and b are small. This double peaking can be explained (see [5]) 
by the multimodality of the Markov binomial distribution, i.e., {/n(^)}fc=o bimodal 
or even trimodal when both a and b are small (see [1 ). 



4. A SYMMETRY PROPERTY OF ASYMMETRIC RANDOM WALK 

The first column of Figure [l] suggests that when a + = 1 the conditional variance is 
symmetric. This is the content of the next proposition. 

Proposition 4.1. For n> 1 and < x < n, if v ~ n and a + b = 1, we have 

Var (5y(n) \ Sxin)=n + x) = Var {Syin) \ Sx{n)=n-x). 

Proof. For ly — tt and a + b = 1, the Markov chain {Zk} is trivial in the sense that {Zk} 
is independent identically distributed. Then the behavior of the particle is a random walk 
starting at position (0, 0) with the following transition probabilities for fc = 0, 1, • • • , n — 1, 

P {S{k + 1) = S{k) + (2, 0)) = ba, P {S{k + 1) = S{k) + (0, 0)) ^ba + a, 
P {Sik + 1) = S{k) + (1, 1)) = P {S{k + 1) = Sik) + (1, -1)) = &/3. 

For convenience we shift the random walk steps by adding (—1,0) to them and let the particle 
start at position (n,0). Then the corresponding transition probabilities of the random walk 
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are given by 

P (^S{k + 1) = S{k) + (1, 0)) = ba, P (^S{k + 1) = S{k) + (-1, 0)) = ba + , 
(^(fc + 1) = Sik) + (0, 1)) = P (^(fc + 1) = Sik) + (0,-1)) = 6/3. 



P 



The result now follows from the following theorem with lo =^ ba, e = ba+a and 7 = (5 = 6/3. □ 

Theorem 4.1. Let (Sn) — iiSn \ Sn'')) be a simple random walk in the plane, i.e., Sq = 
(0, 0) and (Sn) has increments Xk with 

P(A'fc = (1,0)) P(A'fc = (-1,0)) P(A'fc = (0,1)) =7, P(A'fe = (0,-1)) =<5, 

where uj + £ + ^ + 5= 1. Then for a; = 0, 1, 2, . . . , n and — n < y < n, 

P (^(^) = V I ^(1) = x) = P = y I ^(1) = -x) . 

Proof. A similar combinatorial analysis as in the proof of Proposition |3.1| yields that for 
a; = 0, 1, • • • , ri and —n < y < n 

n n -y)- 1^ J-! (j- + :e)! {{n -X- y)/2 - j)\ {{n - x + y)l2 - j')! ' 



_ _ . ^ ^ = ^,^]+x^(n-x-2,+y)/25{n-x-23-y)/2 



^.^0 j! (j + x)\ ((n - X - y)/2 - j)\ {{n -x + y)/2 - j)\ 



iiy = n-x mod 2, and P [sl^^ = a;, S*!^' = y) = P [s^n^ -a;, 5^^^ = y) = if y ^ n - a; 
mod 2. Similarly, 

i! (i + a;)! (n — a; — 2i)! 
7 1 



i=0 



= ^ ^,.V+-(7 + <5)»-2.-._ 

V / ^-^ i! (i + a; ! m, — a; — 2z )! 



^ i! (i + a;)! (n - a; - 2z)! 
Then it is easy to check that for x = 0, 1, • • • , n and — n < y < n 

P (5(1) = X, = y) P (^(1) = -x) = P (^(1) = -X, 5(2) = y) P 
which finishes the proof of the theorem. □ 

5. Forty five degree dispersion 

To fix the strange behavior observed when the dispersion is 2D simple random walk, we 
consider now a well known variant with increments 

P((Xfe,rfc) = (l,j)) =a, P((Xfe,rfc) = (-l,j)) =/3 forj-il, 

where a,/3 > and a + /3 = 1/2. Let [S{n)) be the corresponding kinetic transport process. 
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Proposition 5.1. For n > I and < x < n, —n < y < n, the distribution of S{n) is given 
by 

and the conditional variance is given by 

EL.fc/n(fc)C)(2«)^(2/3)'=-^ 



EL./"(fc)©(2a)-(2/3)'=- 



Var(5Y(n) | Sxin) = 2x) 

Proof. The proof of the proposition is quite similar to that for Proposition |3.1| For < a; < n 
and —n < y < n, we consider the probabihty of the particle being at position (2x,y) at time 
n by conditioning on Kn — k. In order to achieve Sxin) = k + Xi + • • • + Xk = 2x the 
particle needs to move x steps to the right and k — x steps to the left. Also to obtain 
Syin) = Yi + ■ ■ - Yk ~ y the particle needs to move (fc + y)/2 steps up and (k — y)/2 steps 
down. Let j be the number of steps of type (1,1), then the particle moves x — j steps of type 
(1,-1), {k + y)/2 — j steps of type (—1,1), and the remaining k — [j + x — j + {k + y) /2 — j) = 
(k—y) /2— (x—j) steps of type (—1, —1). To ensure the nonnegativity of all of these quantities, 
we have x < k < n and max{0, (y — k)/2 + x} < j < min{x, {k + y)/2}. We conclude by 
using Vandermonde's identity, obtaining that for < x < n, ~n < y < n and k = y mod 2 

P (5(n) = (2., ,) I ir„ = fc) = ''g'' Q Q (^-_^^.) a^^^^^^ 

For the conditional variance Var(5y(n) | Sx{n) — 2x), we still need the probability of the 
particle being at position 2x at time n. In a similar way as above it follows that 

P (SxH = 2x I Kr, - fc) - (2a)-(2/3)'=^^ 

By symmetry E[5y(n) | (Sx(^) ~ 2x] = 0. So 

n 

Var (5y(n) | Sx{n) = 2x) = ^ y^V (5y(n) ^y\Sx{n)^ 2x) , 
and the proposition follows by the following calculation 



n n / , 

y'V{S{n)^{2x,y))^ y' ')(^)"'/^'~' 



y — — n k = y mod 2 

a:;<A;<n 



2 



k=x ^ ^ y = fc mod 2 V 2 / fc^a; V / s^o ^ ^ 

— n<y<n 

Y^f^n(k)[\2arm^-\ 



□ 



Example 5.1. Let n = 50, a ^ jS — 1/4, and the initial distribution v — tt. See Figure]^ 
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Figure 2. The probability mass function of 5(50), a list contour plot of 
its mass function and the conditional variance Var (Sy (50) | Sx (50) = 2x) 
with h' ^ TT,a ^ /3 — 1/4. In the first column a — 0.1, b — 0.9; in the second 
column a = 6 = 0.1; in the last column a ^ b — 0.01. 



Figure [2] suggests that the conditional variance Var(5y(ri) | Sx{n) — 2x) is increasing in 
X. This is indeed the case, and can be shown by using Proposition |5.1| 

Proposition 5.2. The conditional variance Var (5y(n) | Sx{n) = 2x) is an increasing func- 
tion in X. 



Proof. Let g{x) EIU ^/"(fc) © (2a)"(2/3)'-" and h{x) EL. © (2«)"(2/3)'=-^ 
According to Proposition 5.1 Var(5i'(n) | Sx{n) — 2x) = g{x)/h{x). For Q < x < n, 
note that Var (5y(n) | Sx{n) = 2{x + 1)) > Var (5y(n) | Sx{n) — 2x) is equivalent to 5(2: + 
l)h{x) > g{x)h{x + 1), since h{x) > 0. Consider the difference 



g{x + l)h{x) - g{x)h{x + 1) 



k—x-{-l i—x 
n 

> E (^-^) 

k,i=x+l 



X + 1 J \X 



fn{k)fnma/(3f^+\2(3)'+' 



X + 1/ \x 
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where w{k,i) = /„(/c)/„(£)(a//3)2^+i(2/3)'=+^ > 0. Using the symmetry w{k,£) = w{i,k), the 
proposition then follows from 

gix + l)h{x) - g{x)h{x + 1) 



n 



k,e=x+i 



k 

x + 1 



- w{kj)ik~£) 

k,l=x+l 



n 

k,e=x+i 

x + l)\x) U + 1 



x + l 



{k~£f >0. 



w{i, k) 



□ 



Proposition 5.2 seems to suggest that the strange phenomenon in the conditional variance 
(not increasing as in Figure [T]) will disappear if we consider independent increments {X^, Yk) 
where Xk is also independent of Yfc for fc = 1, 2, • • • , Kn- However, we will give a mean zero 
nearest neighbor random walk example, where this phenomenon is still present. 

Let {Xk, Yfc), fc > 1 be independent identical distributed random vectors with distribution 

P {{Xk.Yk) = (1, ±1)) = P ((Xfe, Ffc) = (-1, ±1)) = P {{Xk,Yk) = (±1, 0)) = ^ - 2^ 

where < ^ < 1/4. Obviously, X^ is independent of Yfc. Let S{n) = {S x {n) , Sy (n)) be 
the position of the particle at time n. One can compute its probability mass function and 
conditional variance Var | Sx{n) — x). 

Example 5.2. Let n = 25,a = b = 0.01,^ = 0.2, and v = tt. See Figure\^ 





Figure 3. The probability mass function of 5(25), a list contour plot of 
its mass function and the conditional variance Var (5y(25) | Sx{'2-^) = 2a;) 
with = TT, a = 6 = 0.01, $ = 0.2. 



6. Gaussian dispersion 

It is natural (and common practice in the engineering literature, see e.g., [3]) to consider 
a Gaussian dispersion. Let S(7i) be the position of the particle at time n as described in 
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Section [2] except that (Xk,Yk) now has a Gaussian distribution with independent marginals 

Xk = M{0, 2a), Yfc = 7V(0, 2/3) for a, /3 > 0. 

Proposition 6.1. For n > 1, the distribution /i„ of the random variable S{n) can be written 
as 

= /ri(0)^(o,0) + (1 - /n(0))A«, 

where (5(o.o) is Dirac measure at the point (0, 0) and fin is the distribution of a continuous 
random variable having density function 

/s(„)(.:, .y) = 1^ exp ( - - — j . 

A; — 1 

Moreover, the conditional variance is given by 

Var(Sy(n) | S^N = x) = 2/3(1 - /„(0))— — 

L,fe=iexp(- [x- ky/{4ka))fn{k)/Vk 

Proof. As in [2] one easily finds the characteristic function fs{n) of S(n), and as in [2] one 
shows that (0,0) is the unique atom of /i„, with mass /n(0). Then the density function of 
the continuous part of (!„ follows by inverse Fourier transformation: 

kn)i^^y) = YZJj0)J^ J J {^Sin){ui,U2)^fn{0))e-^<-'^^^+-^y^du,du2 

^ 1-] (0) E/"W(^ / ye-™i-+'=(»i-"''i)e-™^^-'='^"^d7.idu: 



2 



(x — k)^ 

/„(0) ^^■"""^'47rfcV^^"'' V 4ka 4fc/3 



Using symmetry, the conditional variance is given by 

Var(Sy(n) | Sxin) ^ x) ^ E {Si-{n) \ S^(n) = x) = (1 - /„(0)) T ^2 /s(»)(^, y) 

"'-00 fSx(n)(x) 

where fsxin){x) = fs{n){x, y)'^y- The proof is then finished by a straightforward calcu- 
lation. □ 

Example 6.1. Let n = 50, a = /3 = 1/4, anrf i/ie initial distribution v = t: . See Figure^ 



Figure [4] suggests that the conditional variance Var(Sy(7i) | Sx{n) = a;) is increasing in 
x. This is indeed the case, and is the content of the following proposition. 

Proposition 6.2. The conditional variance Var(Sy(n) | Sx{n) = x) is an increasing func- 
tion in X. 
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Figure 4. The density /s(5o), a contour plot of the density and the condi- 
tional variance Var (Sy(50) | Sx(50) = x) with ly ^ n, a ~ f3 = 1/4. In the 
first column a = 0.1, b = 0.9; in the second column a = b = 0.1; in the last 
column a — b — 0.01. 



Proof. Since by Proposition 6.1 Var (Sy(n) | Sx{n) = a;) is a differentiable function, we only 
need to show that its first derivative is nonnegative. For simplicity, we write for 



9{x) = ^/„(fc)v^e" 



k = l 



So Var(Sy(n) | Sx(n) = a;)' = ^^^jf^{g' {x)h{x) - g{x)h'{x)). Note that 



g'{x)h{x) - g{x)h'{x) 
1 ^ /„(fc)/„W 



2a ^ \/fc 



1 



e — , — (x-k) 
k 



2a 



k,e=i 
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where w{k,£) = -^^^ ^^^^e" ' - ' > 0. Using the symmetry w{k,£) = w{£,k), the 
proposition then follows from 
g'{x)h{x) -g{x)h'{x) 

= ^[ E w{k,£)^-^{x-k)+ w{i,k)'^{x-e) 

1 " ^ i f i.^2 

= e; E - ml - 1) - (i - 1)1 = I; E "ftfl^ > 0. 



□ 
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